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Multi-boson symmetrization effects on two-particle Bose-Einstein interferometry 
are studied for ensembles with arbitrary multiplicity distributions. In the general 
case one finds interesting residual correlations which require a modified framework 
for extracting information on the source geometry from two-particle correlation 
measurements. In sources with high phase-space densities, multi-boson effects 
modify the Hanbury Brown-Twiss (HBT) radius parameters and simultaneously 
generate strong residual correlations. We clarify their effect on the correlation 
strength (intercept parameter) and thus explain a variety of previously reported 
puzzling multi-boson symmetrization phenomena. For event ensembles of (ap- 
proximately) fixed multiplicity, the residual correlations lead to a minimum in the 
correlation function at non-zero relative momentum, which can be practically ex- 
ploited to search, in a model-independent way, for multi-boson symmetrization 
effects in high-energy heavy-ion experiments. 

1 Introduction 

Heavy-ion collisions at high energies produce large numbers of secondary par- 
ticles. Even if these particles are created independently ("chaotic source") and 
are free from final state interactions, the measured n-body momentum spectra 
still do not simply factor into a product of single-particle spectra. Instead, 
symmetrization of the many-body wave function with respect to the exchange 
of identical particles leads (in the case of bosons) to so-called Bose-Einstein 
correlations. 

Bose-Einstein correlations affect the measured cross sections whenever 
more than one boson occupies an elementary phase-space cell. This is the 
basis of Bose-Einstein or intensity interferometry: for a source of size R, sym- 
metrization effects lead to an enhancement of, for example, the two-particle 
cross-section at relative momenta q < 1/R. The theory which allows to extract 
the source size R from two-jparticle correlation data assumes, in its simplest 
and usually employed form,EJ that the effect is dominated by symmetrization 
with respect to an exchange of the two observed particles, and that the ad- 
ditional symmetrization with respect to all the other identical bosons in the 
multiparticle final state can be neglected. 

This is only true if the phase-space density of the source at freeze-out 
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(i.e. at point of last interaction) is sufficiently small. In regions of high phase- 
space density, multi-boson symmetrization effects change the shape even of 
the single-particle spectra and two-particle correlations. They lead to a clus- 
tering a low momenta in momentum space and at low separation in coordinate 
space. For example, in an infinite medium they turn an exponential momen- 
tum distribution of distinguishable elementary pion production sources into a 
Bose-Einstein distribution for the observed pion momenta. B'cl And in two-pion 
interferometry, the effective source extension extracted from the two-pion cor- 
relation function is smaller thaii-the r.m.s. width of the spatial distribution 
of the individual pion emitters.QB 

In addition to these generally accepted phenomena, multi-boson sym- 
metrization effects also cause, other, somewhat more confusing effects. For 
example, it has been notedo^Q that in ensembles with fixed total pion mul- 
tiplicity multi-pion effects reduce the intercept of the two-pion correlation 
function at vanishing relative momentum. This has been interpretedO as a 
sign of "effective coherence" in the source which becomes significajit if its 
phase-space density is large and which will lead to a "pion laser"!] once a 
critical phase-space density is exceeded. This interpretation will be shown 
to be incorrect: The intercept reduction happens even for completely chaotic 
sources without any phase coherence, and even the formaition of a Bose con- 
densate is not accompanied by "pion lasing" . We foundQ that the reduced 
intercept reflects residual correlations among the pions which are not directly 
related to the Hanbury Brown-Twiss correlations which are exploited in Bose- 
Einstein interferometry. These residual correlations are generic in the sense 
that they vanish only for very specific cases where the multiplicity distribution 
of the ensemble takes on a particular form.B When present, they render the 
extraction of the source size from Bose-Einstein interferometry more difficult. 

2 Multi-boson formalism 

Due to space reasons, I can give here only a very rudimentary account of 
the mathematics; alLteplmical details can be found elsewhere.B Following ear- 
lier work by others ,EH3'LrEij the source is parametrized by a density operator 
constructed out of a superposition of n-particle states emitted from classical 
source currents which themselves are Gaussian wave packets. The number of 
these source current wave packets, their centers in phase-space, and the pion 
multiplicity n all have completely arbitrary (although normalizable) distri- 
butions. The source currents emit pions with random phases which ensures 
complete source chaoticity. For techiiijqai_reasons, the n-pion states used in 
the construction are not normalized,l3i3'ElEil but the density operator is. Cal- 
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culating this normalization is the only part of the problem which cannot be 
done analytically even for the simple case of nonrelativistic, instantaneous 
Gaussian—wave packets, but requires the numerical evaluation of a recursion 
relation .113 |^ 

As found by PrattiZl, all n-particle cross sections can be expressed in terms 
of the so-called "ring integrals" 

Gi{Pi,P2) = J Gi{Pi,ki)dki Gi(fci,fc2) • ■ ■*i_iGi(fci_i,p2) , (1) 

where dk = (Pk/ E^. They describe the effect on the two-particle exchange 
amplitude (with respect to the interchange of momenta pi and P2) of 
the symmetrization with respect to 1 additional particles with momenta 
They can be generated recursively from the elementary two- 
particle exchange amplitude Gi{pi,p2), the Fourier transform of the source 
Wigner density g(x, K): 

Gi(pi,P2)= [d''xg{x,'^{p^+p,))e-^^P^-P^y-^ f d^x g [x, K) e~^'i \ (2) 



For a non-relativistic spherical Gaussian Wigner density with spatial width R 
and momentum spread A, which emits paridjdes instantaneously '-^ 5{t), the 
ring integrals can be calculated analyticallyflB 



GrM.P.) - exp (^-^ - ' 

Rl = anR\ Al^anA\ v = 2RA > 1 , (4) 

One sees that the higher order ring integrals are Gaussians with reduced 
widths Rn and A„ in both coordinate and momentum space, reflecting the 
bosonic clustering at low momenta and relative distances mentioned in the 
Introduction. The strength of these clustering effects is controlled by the 
single parameter v = 2RA which measures the phase-space volume of the 
source. The above equations remain correct for expanding Gaussian sources 
(after suitable redefinitionil of the parameters R and A), and they are easi- 
ly generalized for Gaussians with different widths along the three Cartesian 
directions. 
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3 Bose-Einstein and residual corelations 

For this very general class of models, the two-particle correlation function can 
be written down explicitly as follows: 

Er^=i hi+j Gi{pi,pi) Gj{p2,P2) 



C2{Pl,P2) 



J2Z = 1 ^« ^3 '^i ' ^1 ) '^0 , P2 ) 

hi+j Gi {pi , P2) Gj (p2 ,piy 



1 



T,Z = 1 ^^+'J Gi{Pl-,Pl) Gj{p2,P2) J 

^Gt'\q,K)(l + R2{q,K)). (6) 

The sums over i and j take into account symmetrizations with respect to an 
increasing number of other particles in the final state. The weights hi are 
related to the measured pion multiplicity distribution p„ by 



h^ = Y,Pn ^'j^C , (7) 



Lo{n — i) 
uj{n) 



where uj{n) is (up to trivial factors) the factor required to normalize the n- 
particle density operatorH and must be obtained by numerically solving the 
recursion relationlla 

1 " 

cj(n) = — > CiUj{n — i), (8) 

i=l 

with Lu(0) = uj{l) = 1 and Gi = J dpGi{p,p) ~ c„an^^. Obviously, the Lu(n) 
and the weights hi depend non-trivially on the phase-space volume factor v. 

For dilute systems {v 3> 1), the higher order ring integrals can be ne- 
glected, and all the sums in are dominated by their fist terms. In this case 
the prefactor Cj'^"'''' in reduces to a nornialization constant J\f = h2/hf, 
and the correlator assumes its standard form:EJ 

C2[pi.,P2) [1+ ^N^/^ ■ (9) 

In general, however, the prefactor is a non-trivial function of both the pair 
momentum K and the relative momentum q. These residual correlations do 
not involve the exchange of momenta pi and p2 and thus are not related 
to Bose-Einstein interferometry; they cannot be used to extract the source 
radius R. They do, however, modulate the correlation function and modify 
its intercept at g = 0. A typical example is shown in Figure 0. One sees 
the genuine (Gaussian) Bose-Einstein correlation function R2{K, q) sitting on 
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Figure 1. Two-particle correlation function, indicating separately the Bose-Einstein correla- 
tor i?2 and the residual correlations Cj'^^'^' , for a fixed pion multiplicity n = 10, a phase-space 
volume factor v = 2RA = 2.56, and vanishing reduced pair momentum k = K/ A = 0. 



top of the residual correlation; the latter starts out at the value Af = h2/h\ at 
q = oo but decreases to a much smaller value at q = 0. 

The source size signal is contained in the function R2 in (||). Note that 
i?2 goes to 1 as q ^ 0; thus at g = the exchange term is as large as the direct 
term, implying a fully chaotic source (as it should be by construction). The 
reduced intercept at q — thus has nothing to do with coherence. 

The residual correlations can be positive or negative, enhancing or de- 
creasing the intercept at q = above or below the value 2J\f. Negative residual 
correlations appear, for example, in ensembles with fixed pion multiplicity or 
with a Poissonian multiplicity distribution; on the other hand, a Bose-Einstein 
multipliciy distribution can lead to large positive residual correlations.u The 
residual correlations disappear in the dilute limit, v —^ oo, but generically 
become strong for large phase-space densities or small values of w._At the mo- 
ment there are only two special multiplicity distributionS|_knowra for which 
residual correlations are completely absent, for all allowedllj values of v. 

The limit u — > 1 is particularly interesting. In this limit the source is 
as small as allowed by the uncertainty principle. The result are vanishing 
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HBT radiija and the correlation function becomes completely flat (i.e. q- 
independent) . [Note that in this limit also the residual correlation function 
6*2'^°''^ become g- independent.!] Superficialhj-his looks like the correlation func- 
tion from a coherent state ("pion Iaser"l3l3) which shows no Bose-Einstein 
correlations. This impression is, however, quite misleading: The correlation 
function still has the structure (j^) , with an exchange term R2 = I which is as 
large as the direct term (i.e. the source is still fully chaotic). Note that this 
limit can be taken while keeping the pion multiplicity distribution fixed, i.e. 
the latter does not automatically become a Bose-Einstein distribution when 

4 Epilogue 

Does all this matter for our real lives as heavy-ion physicists? Most likely 
not. An analysis of average pion phase-space, densities at freeze-out from 
heavy-ion collisions at AGS and SPS energiesEj yields universally low values 
for all considered collision systems and beam energies. First results from 
RHIC appear tOj-cpnfirm this even at the much higher collision energy of 
y/s — 130AGeVJl3 Apparently pions stop interacting (freeze out) only after 
the phase-space density has dropped below a critical value which is low enough 
for multi-boson symmetrization effects to remain unimportant. This is good 
news: it means that we can continue to use the standard HBT formalismld to 
extract source sizes from two-pion correlation data. It would be good, though, 
to confirm this by directly checking for (the absence of) residual correlations 
in a high-statistics analysis of two-particle correlation data from heavy-ion 
collisions. This can be done by selecting heavy-ion collisions producing an 
(approximately) fixed number of, say, tt"*" mesons and searching for a nainimum 
at non-zero q in the (Coulomb corrected) 7r+7r+ correlation function.El At the 
same time, it may be interesting to explore other fields of application for the 
elegant multi-boson symmetrization techniques presented here. 
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